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FMOS near-IR spectroscopy of Herschel selected galaxies: 
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ABSTRACT 

We investigate the properties (e.g. star formation rate, dust attentuation, stellar mass 
and metallicity) of a sample of infrared luminous galaxies at z ~ 1 via near-IR spec- 
troscopy with Subaru-FMOS. Our sample consists of Herschel SPIRE and Spitzer 
MIPS selected sources in the COSMOS field with photometric redshifts in the range 
0.7 < z phot < 1.8, which have been targeted in 2 pointings (0.5 sq. deg.) with FMOS. 
We find a modest success rate for emission line detections, with candidate Ha emission 
lines detected for 57 of 168 SPIRE sources (34 per cent). By stacking the near-IR spec- 
tra we directly measure the mean Balmer decrement for the Ha and H/3 lines, finding 
a value of (E(B-V)) = 0.51±0.27for (L m ) = 10 12 L Q sources at (z) = 1.36. By com- 
paring star formation rates estimated from the IR and from the dust uncorrected Ha 
line we find a strong relationship between dust attenuation and star formation rate. 
This relation is broadly consistent with that previously seen in star-forming galaxies 
at z ~ 0.1. Finally, we investigate the metallicity via the N2 ratio, finding that z ~ 1 
IR-selected sources are indistinguishable from the local mass-metallicity relation. We 
also find a strong correlation between dust attentuation and metallicity, with the most 
metal-rich IR-sources experiencing the largest levels of dust attenuation. 

Key words: galaxies: evolution, submillimetre: galaxies 
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1 INTRODUCTION 

Accurate measurements of the characteristic properties of 
galaxies; star formation rate (SFR), metallicity and stellar 
mass, are central to our understanding of their evolution. 
Significant progess in our ability to measure these proper- 
ties in distant (z > 1) galaxies has been made in the last 
two decades. Deep surveys with the Hubble space telescope 
have enabled the star formation rates of large numbers of 
galaxies up to z ~ 7 to be measured (e.g. Madau et al. 
1996; Bunker et al. 2004; Bouwens et al. 2006; Bouwens et 
al. 2009., McLure et al. 2010). Large scale optical and near- 
IR spectroscopic surveys have targetted key emission lines 
allowing the study of gas metallicity to z < 4 (e.g. Tremonti 
et al. 2004; Erb et al. 2006; Mannucci et al. 2009; Zahid et al. 
2011; Cresci et al. 2012). Finally, deep optical and near-IR 
photometric surveys have allowed an accurate assessment of 
the stellar mass contained in galaxies, and its build-up with 
redshift out to z ~ 5 (e.g. Fontana et al. 2006; Ilbert et al. 
2010; Caputi et al. 2011). 

While these advances have re-shaped our understanding 
of galaxy formation and evolution, they typically rely obser- 
vations in a single wavelength window i.e. UV, optical/near- 
IR, far-IR, etc. Meanwhile at low and intermediate redshifts 
it is becoming clear that large scale, multi-wavelength, stud- 
ies of galaxies are needed to determine unbiased estimates of 
their properties. In the case of SFR estimates, where the cor- 
rections for dust attenuation tend to be large, comparisons 
of UV and IR SFR estimates at z ~ (Hao et al. 2011), 
z ~ 1 (Buat et al. 2010) and z ~ 2 (Reddy et al. 2011) 
show that widely used in band (i.e. in the same waveband 
as the SFR estimate) dust attenuation estimators (e.g. the 
UV continuum slope; Meurer et al. 1999) have large errors 
(<51og 10 SFR ~ 0.3 dex) and can have significant systemic 
biases for certain populations of galaxies (low SFR spirals, 
ULIRGs lO 12 Z/0 and young starbursts). By comparison SFR 
estimators based on combinations of the far-IR, Ha line or 
radio have much smaller errors (8 log 10 SFR ~ 0.1 dex; Ken- 
nicutt et al. 2009; Hao et al. 2011) and are universally valid. 

Measurements of other galaxy properties also bene- 
fit from a multi-wavelength approach. Metallicity estimates 
from single tracers e.g. the N2 or R32 methods (Pettini & 
Pagel 2004) can disagree by up to A[log 10 (O/ff)] = 0.7 dex 
(Kewley & Ellison 2008) . Stellar mass estimates obtained via 
the fitting the stellar population models require multi-band 
observations in the optical and near-IR to be reliable; omit- 
ting near-IR observations introduces an error of ~ O.ldex to 
stellar mass estimates at 2 ~ 1 (Pozzetti et al. 2007; Ilbert 
et al. 2010). 

In order to put galaxy evolution at high— z on a firm 
footing multi-wavelength observations at the same rest- 
frame wavelengths as our low-z benchmarks (e.g. SDSS, 
Spttzer, IRAS) for a large number of high-z galaxies is 
needed. This necessitates wide-field imaging and spectro- 
copy in the IR. 

Here we investigate the rest-frame optical-to-far IR 
properties of a sample of Hersche^\ (Pilbratt et al. 2010) 

Herschel is an ESA space observatory with science instruments 
provided by Principal Investigator consortia. It is open for pro- 
posals for observing time from the worldwide astronomical com- 
munity. 



sources which were targetted for near-IR spectroscopy with 
FMOS (Kimura et al. 2010). The key goal of this work is 
to determine the key galaxy properties (SFR, dust atten- 
uation, stellar mass and metallicity) between a sample of 
high— z (0.8 < z < 1.7), IR luminous (> 10 11 L Q ) sources 
using the same tracers commonly used for low-2 samples. 
In this way we can be sure that our results are fully consis- 
tent (in terms of both calibration and selection effects) with 
those at low-z. The datasets used in this work are described 
in Sj2] £13.11 presents the detection rate of Ha, ^13.21 the ag- 
gregate near-IR spectral properties and £13.31 a comparison 
of the star formation rates from the IR and Ha line. In > )3.4I 
we investigate the stellar mass and metallicity of our sample 
and, finally, Sj4] summarises our conclusions. Throughout we 
assume a ACDM cosmology with J1a = 0.7, fi m = 0.3 and 
H = 70kms- 1 Mpc- 1 . 



2 DATA 

2.1 Pre-existing COSMOS data 

The starting point for this work is the SPIRE observations 
of the COSMOS field (Scoville et al. 2007) taken as part of 
the Herschel Multi-tiered Extragalactic Survey (HerMES; 
Oliver et al. 2012). The SPIRE instrument, its in-orbit per- 
formance and its scientific capabilities are described by Grif- 
fin et al. (2010); its calibration methods and accuracy are 
outlined in Swinyard et al. (2010). Here we make use of 
SPIRE maps as described in Levenson et al. (2010). At 
the time of writing HerMES observations of COSMOS cover 
~ 4.8 sq. deg. to a la instrumental noise of ~2mJy/Beam 
at the three SPIRE wavelengths of 250, 350 and 500 /mi. 

As the SPIRE data offer an instrumental noise signifi- 
cantly lower than the confusion noise (~ 6mJy; Nyugen et 
al. 2010) we make use of prior source positions from higher 
angular resolution data to extract SPIRE photometry. The 
MIPS 24 /im channel is the most obvious prior for SPIRE 
data as it offers a significant improvement in angular res- 
olution (6 arcsec FWHM for MIPS 24 /im vs. 18.6 arcsec 
FWHM for SPIRE 250 /im) while also being able to account 
for > 80 per cent of the flux at SPIRE wavelengths at the 
24 /im depths now available in a large fraction of HerMES 
fields (Bethermin et al. 2012; Oliver et al. 2012). 

To construct our prior catalogue for SPIRE photome- 
try we begin with the MIPS 24 /im imaging from the Spitzer 
COSMOS survey (Le Floc'h et al. 2009). Here we make use 
of the publicly available imaging, performing source extrac- 
tion via the STARFINDER IDL package (Diolaiti et al. 2000). 
The resulting catalogue covers ~ 2.1 sq. deg. and has a typ- 
ical lcr sensitivity of a — 15 /iJy. 

In order to provide the most accurate positional infor- 
mation, for both our SPIRE photometry and FMOS fibre 
positioning, we cross-match our 24 /im catalogue to the pub- 
lically available HST ACS /F8i4W-band catalogue of Leau- 
thaud et al. (2007). This catalogue covers 1.64 sq. deg. to 
a limiting magnitude of /fsi4w < 26.5. In addition to im- 
proving the positional accuracy this matching helps elim- 
inate spurious 24 /im sources produced by artifacts in the 
image, in particular those located close to bright (> lmjy) 
sources. Of the 35,914 24 /im sources located well within the 
ACS i-band coverage, 33,071 (92 per cent) have /F8i4w-band 
counterparts within 2 arcsec. 
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SPIRE photometry is performed using the Jp8i4w-band 
positions of the 24 /xm sources as a prior, following the algo- 
rithms described in Roseboom et al. (2010) and Roseboom 
et al. (2012). All > 3a (~ 60 /iJy) 24/im sources are con- 
sidered as potential SPIRE counterparts. Using the residual 
map statistics we estimate that our prior-driven SPIRE cat- 
alogue reaches a typical point source sensitivity of crtot =2.7, 
3.5, and 3.2 mjy at 250, 350 and 500 /im, including the con- 
tribution from source confusion. 

To complete our multi-wavelength COSMOS dataset 
we add multi-band optical/near-IR data and photometric 
redshifts from the catalogue of Ilbert et al. (2009). This 
catalogue is limited to < 25, and hence our HerMES- 
COSMOS sample is similarly restricted. 

2.2 FMOS observations and emission line 
measurements 

IR-selected sources were targeted in two pointings (0.5 
sq. deg.) located within the COSMOS field with FMOS as 
part of the GTO program. The FMOS instrument (Kimura 
et al. 2008) consists of 400 1.2-arcsec diameter fibres which 
can be placed within a 30-arcmin diameter field of view. We 
used the low-resolution mode (R ~ 600) , allowing instanta- 
neous coverage of both the J and H band (0.9 < A < 1.8 fiva), 
with cross-beam switching i.e. two fibres for each target; one 
placed on the sky and one on the target, with the target/sky 
"switched" between them at regular intervals. 

Potential targets for FMOS fibre allocation were se- 
lected from our HerMES-COSMOS parent catologue by re- 
quiring a photometric redshift in the range 0.65 < z p hot < 
1.75 from the catalogue of Ilbert et al. (2009). This restric- 
tion was introduced to ensure that the Ha line was likely 
within the FMOS wavelength coverage. 

Fibre allocation preference was given to sources de- 
tected at both 24 (jm and 250 /im (> 3<7 to t), followed by 
24 /im only sources. As well as science targets, a number 
(typically 2-4) of 2MASS (Skrutskie et al. 2006) selected 
stars were included in the observations for flux calibration 
purposes. 

The first of our FMOS pointings was dedicated to solely 
HerMES-COSMOS targets, while for the second pointing 
(2010 November 24 and 25) we shared fibres with the 
evolSMURF project (Bunker et al., in prep). While the split 
between the samples was roughly 50-50, this was aided by 
the overlap between the samples (42 sources). 

Table [1] details the exposure times, together with the 
number of 24 and 250 /im detected (henceforth refered to as 
24 n 250 /im), and 24 /im only sources in each pointing. In 
total 241 fibres were allocated to IR-selected sources, with 
four sources appearing in both configurations, resulting in 
237 unique targets (168 unique 24R250 /im targets). All data 
were reduced using the standard FMOS pipeline (Iwamuro 
et al., 2011). 

Emission lines were identified in the 2D reduced frames, 
after flux calibration, via a semi-automated procedure. At 
each pixel the line profile was fitted to the surrounding 9x9 
pixels. The pixel scale was 5A in the spectral direction and 
0.13 arcsec in the fibre direction. We only considered pix- 
els in the wavelength ranges 1.1-1.36 /im and 1.42-1.7/im. 
Pixels within 5 A of an OH line were excluded from consid- 
eration. The line profile was assumed to be Gaussian with 



Table 1. Summary of FMOS observations 





Date 


T 

cxp 


-^24 firn 


^24n250 iim 


PI 


2010 Nov. 22 


8 x 900 s 


136 


102 


P2 


2010 Nov. 24 & 25 


14 x 900 s 


105 


67 




Total (Unique) 




241 (237) 


169 (168) 



FWHM= A/600 A in the spectral direction, and 6.9 pixels 
in the fibre direction. The noise was estimated by taking 
the variance of all illuminated pixels at that wavelength on 
the detector. Regions where the noise is exceptionally high 
(> 10/xJy pixeP 1 ) were excluded. The local continuum was 
estimated by taking the median pixel value in a window of 
20 pixels (200 A), excluding the closest 7 pixels. For each fit 
the line signal-to- noise ratio (SNR), peak SNR and correla- 
tion coefficient between the line profile and 2D spectrum was 
measured. The line SNR (crn nc ) was calculated via Eqn.[T] 

-, -1/2 



<TU, 



=5> 



(1) 



where di is the pixel intensity at position i, Ci is the contin- 
uum at pixel i, Pi the line profile at position i, and cr; the 
noise estimate at position i. 

The peak SNR is defined as the ratio of the peak flux 
density, taken to be the mean flux density in a 3x3 pixel 
window less the local continuum, to the standard deviation 
of the surrounding pixels in the spectrum. 

Finally the correlation coefficient (pn n o) is calculated 
Eqn.H 

Pline = (4) 

where op and Od are the standard deviation of the line pro- 
file and data values, respectively. 

Line fits which have (Jii no > 4, o" pca k > 2.5 and pii nc > 20 
were considered as candidate emission lines. 

All candidate emission lines that have a wavelength 
within the range (1 + z phot - 0.16) x 6563.4A< A < (1 + 
zphot + 0.16) x 6563.4A were considered to be Ha. The win- 
dow of 8z — 0.16 equates to 4<7 p hot-z for the typical photo-z 
error quoted by Ilbert et al. (2009) at z ~ 1. 

A total of 85 candidate Ha emission lines were found 
from the sample of 237 unique 24 /im targets. We assessed 
the reliability of our line identification technique in two 
ways. Firstly the line identification was repeated, but with 
the proposed redshifts (and hence wavelength search win- 
dow) shifted. To ensure that the mock search windows are 
sufficiently far away from real lines, but still within the wave- 
length coverage of FMOS, sources at z p hot < 1-1 were given 
Zmock = z phot + 0.16 + 8, while those at z phot > 1.1 were set 
to Zmock = z P hot — 0.15 — 8, where 8 is a random number 
between and 0.15. As a result of this process seven lines 
were identified, giving an estimate of the false positive line 
detection rate of 7/85 or 8±3 per cent. No false lines were 
returned with o"ii no > 8. 

We compare the redshifts, as determined by the wave- 
length of our candidate Ha lines, and the known spectro- 
scopic redshifts. From our sample of 85 Ha line emitters, 
28 are found to have reliable (z qua i > 3) spectroscopic red- 
shifts in the zCOSMOS bright (i < 22.5) sample of Lilly 
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Table 2. Completeness of line identification process as a function 
of flux. Completeness is estimated via injection of mock emission 
lines into the 2D spectra at random wavelengths/redshifts. Note 
that this is for flux contained within the 1.2 arcscc diameter fibre 
of FMOS; no aperture effects are considered. 



Line flux Completeness 
PI P2 

10 — 16 ergscm -2 s — 1 . percent percent 



0. 5 4.5 1.4 

1. 23.2 25.3 
1.5 52.9 57.8 

2. 65.8 67.8 
2.5 71.0 72.1 

3. 74.2 74.7 
3.5 77.4 77.6 

4. 77.6 78.7 



et al. (2007). Of these, 27 (96 per cent) are found to be 
within 8z — 0.01 of our assumed Ha redshift. No incorrect 
redshifts are found amongst the 11 sources with candidate 
Ha lines at au nc > 8 and a spectroscopic redshift from Lilly 
et al. (2007). While this result is encouraging, it is likely 
that the reliability of our line identification is a function of 
brightness, and the zCOSMOS bright sample is limited to 
*ab < 22.5; roughly 54 per cent (46/85) of our candidate 
Ha line emitters have 22.5 < i\ B < 25. 

Via a similar process we can estimate the complete- 
ness of our line identification process. Assuming the ran- 
dom redshifts, z moc k, described above, we inject mock emis- 
sion lines into our data at a wavelength corresponding to 
1 + Zmock) 6564. 3A. We estimate the completeness for the 
two pointings independantly. Table [5] details the complete- 
ness (i.e. the ratio of sources detected to those injected) 
as a function of flux. No aperture correction is assumed, in- 
jected line fluxes are considered to be those contained within 
the 1.2 arcsec diameter fibre of FMOS. The completeness 
never reaches 100 per cent as lines at certain wavelengths 
will always be undetectable due to the gap in wavelength 
coverage from 1.36-1.42 /jm, as well as the masking of OH 
sky lines. In total 125 OH lines are suppressed which, com- 
bined with the gap due to atmospheric absorption, remove 
24 per cent of the potential wavelength coverage. It can be 
seen that our identification process reaches this maximum 
level above a line flux of /hq Z 4 x 10 -16 ergscm _2 s _1 , 
while we are ^ 50 per cent complete above a line flux 
of /hq k, 1.5 x 10 -16 ergscm _2 s _1 . The completeness es- 
timates for the second pointing are marginally higher at 
/hq S 3x 10~ 16 ergscm~ 2 s -1 due to the increased exposure 
time (12600 sec for the second pointing vs. 7200 sec for the 
first). Note that these line fluxes consider only the flux con- 
tained within the 1.2 arcsec diameter fibres used on FMOS; 
no aperture correction has been applied yet. 

For each Ha line emitter we attempt to measure the flux 
of the neighbouring [Nil] 6584 line. Line fluxes are estimated 
in the same manner as Ha, but with the central wavelength 
fixed at (1 + z) 6584A. Of the 85 Ha line emitters; 33 also 
have [Nil] at SNR > 3. 

Finally, both Ha and [Nil] line fluxes are corrected for 
the limited aperture and unknown stellar continuum. The 
aperture correction is determined via the ratio of 7fsi4w 



flux within the FMOS fibre (1.2 arcsec diameter) to that 
within the Kron radius. The typical aperture correction is 
~ 2-3. For the continuum correction, the near-IR spectra are 
not deep enough to detect the continuum emission near the 
Ha line. Thus, continuum emission at the wavelength of the 
line is estimated, and removed, using the available broad- 
band optical and near-IR imaging from the COSMOS sur- 
vey (Capak et al. 2007, McCracken et al. 2010). For sources 
where the line lies at A < 1.4 /im we use the J band magni- 
tude, and the [z + — J] colour to estimate the continuum flux. 
For those lines where A > 1.4 /im the K band magnitude and 
[J — K] colour are used. In all cases we assume that the Ha 
line is coincident with stellar absorption of equivalent width 
EW = 4.4 (Moustakas & Kennicutt 2006). 

Appendix [X] contains the line fluxes and photometric 
properties for our 85 Ha line emitters. While we list the 
full catalogue of Ha line detections resulting from our ob- 
servations, in the following analysis we consider only the 57 
Ha line emitters which have robust detections at 250 /im. 
This restriction is implemented as 24 fim alone is not a good 
tracer of total IR luminosity at these redshifts (Elbaz et 
al. 2010). Enforcing 250 /im detections also helps to min- 
imise contamination by AGN activity (Hatziminaoglou et 
al. 2010). 



2.3 SDSS comparison sample 

To establish a low redshift baseline for our z ~ 1 mea- 
surements we identify a sample of Spitzer 160 /im selected 
sources from the SWIRE survey (Lonsdale et al. 2003) in the 
Lockman Hole field with SDSS DR6 (Adelman-McCarthy et 
al. 2008) spectroscopy. Sources are required to be brighter 
than 50 mjy at 160 fim and to lie at z < 0.4. The 160 /im 
photometry is complemented with SWIRE MIPS 24, 70 fjtm 
data, and SPIRE photometry at 250, 350 and 500 /im from 
HerMES. Ha line fluxes are calculated using the GAN- 
DALF package (Sarzi et al. 2006), with aperture correc- 
tions based on the ratio of the SDSS Petrosian-to-fibre r- 
magnitudes. We retain sources with a peak line flux density 
to noise ratio of greater than three. The typical limiting Ha 
line flux (3cr) is 5.8 x 10 16 erg cm -2 s _1 . Our final low-z sam- 
ple consists of 156 160 /im-detected sources with reliable Ha 
line flux measurements, at a mean redshift of z — 0.1. 



3 RESULTS 

3.1 Ha detection rate for IR-selected sources 

The raw Ha detection rate for IR-sources (24R250 fim) is 
57/168 (34±4 per cent). In Fig. [T] we present the Ha detec- 
tion rate as a function of 250 /im and 24 /im flux densities. 
While the Ha detection rate appears to be insensitive to 
250 /J,m flux density, a modest increase in the detection rate 
is seen for sources with S24 M m > 200 fi Jy. 

Fig- EH shows the redshift distribution of IR-sources tar- 
geted with FMOS, and those with Ha line detections. The 
contrast between the two distributions highlights the visi- 
bility of Ha with FMOS as a function of redshift. Prior to 
FMOS observations, only photometric redshifts were avail- 
able for the vast majority of targets. Thus while targets were 
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Figure 1. Ha line detection rate of FMOS targeted sources as a 
function of 24 /an and 250 fim flux density. The detection rate is 
relatively insensitive to 250 (im flux density, while a modest gain 
is seen for sources with 524 pm > 200 fi Jy. 



photo— z selected to be in the redshift range where Ha is vis- 
ible, some fraction of sources were expected to have Ha fall 
at wavelengths outside the FMOS wavelength range. 

Despite this the observed detection rate compares well 
with what would be predicted given the known complete- 
ness f i|2.2[) . For each detected source we calculate the com- 
pleteness at that line flux, ci, by interpolating the values in 
Table [5] If the Ha line fluxes of the undetected sources are 
distributed similarly to the detected ones (i.e. the bulk of 
the incompleteness is due to OH sky lines) then the sum 
y] , 1/c; should be equal to the number of targetted objects. 
For our 24H250 /im sample l/c; = 177; close to the ac- 
tual number of sources targetted (168). 

We can compare our detection rate for IR-selected 
sources with FMOS to optically selected sources with VLT 
VIMOS. For the VVDS-DEEP survey (Le Fevre et al. 2005), 
the detection rate for i-selected sources with similar magni- 



tudes (median for our sample is 



23.5) and redshifts 



[z ~ 1.2) to our sample is ~ 70 per cent (for 4.5 hr expo- 
sures; Ilbert et al. 2005). While our FMOS detection rate 
is almost one-half of this, IR-sources have quite a low areal 
density and cannot make the most of the large multiplex 
of VIMOS. For our parent sample of 24 & 250 detected, 
*ab < 25 and 0.65 < z p hot < 1-75 sources the areal density is 
~ 1500 per sq. deg. Using these numbers as a guide, the ex- 
pected number of redshifts recovered in a single 2 hr FMOS 
pointing (0.36x200 fibres equals 72 sources), is comparable 
to the number expected from a single 4.5 hr VIMOS point- 
ing (0.7x91 targets equals 64). Thus for Herschel selected 
targets FMOS is a competitive facility for redshift recovery 
in the range 0.7 < z < 1.8. 




W/.////////////////A 
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Figure 2. Redshift distribution of FMOS targeted sources and 
sources with Ha line detections. For FMOS targets we use pho- 
tometric redshifts from Ilbert et al. (2009). For detected sources 
the redshift implied by the location of the Ha line is used. 



3.2 Composite spectrum for Ha detected sources 

While our near-IR spectra are of sufficient quality to ro- 
bustly measure fluxes for bright emission lines, very little 
additional information can be extracted from the individual 
spectra. This is partly due to the low SNR, but also be- 
cause of the large number of OH sky lines in the near-IR, 
which significantly limit the wavelength coverage. However, 
using the spectral coverage unaffected by OH sky lines from 
each of our Ha-detected sources we can build a composite 
spectrum across a reasonably wide and continuous range of 
rest-frame wavelengths. 

Fig. [3] shows the composite rest-frame spectrum in the 
region of Ha and H/3 produced from our sample of 57 
sources. The composite spectrum is produced by first adopt- 
ing a grid in rest-frame wavelength and calculating the con- 
tribution of each observed near-IR spectrum via a Gaussian 
kernel. The FWHM of the Gaussian kernel is taken to be the 
spectral resolution, A/600 A. Pixels within 10 A of an OH- 
suppressed line are excluded, as are those which lie outside 
the wavelength range 1.1-1.36 /im or 1.42-1.7 /im. The vari- 
ance in the composite spectrum is estimated via jack-knife 
resampling of the sources contributing to each wavelength 
bin. The absolute flux calibration of the composite is cor- 
rected by comparing the continuum level in the composite 
to the mean broadband photometry in the J and K bands. 

Several well-known spectral lines can be seen in addition 
to Ha: [Nil] 6584, although this is blended with Ha; a blend 
of the [SII] doublet at 6716 and 6731 A; and weak signatures 
of [OIII] 5008 and H/3. 

The line ratio Ha/ [Nil] 6584 compared to the line ra- 
tio [OIII] 5008/H^ is often used as diagnostic of AGN ac- 
tivity (Baldwin, Phillips & Terlevich 1981; BPT). Mea- 
suring the strengths of these lines from Fig. [3] we deter- 
mine /[nii] 6584 //hq = 0.24 ± 0.07, while /[oiii] scW/h.s = 

1.3 ± 0.6. Comparing these ratios with the BPT diagnostic 
plots suggests that, on average, the primary origin of emis- 
sion lines, and by proxy the IR-luminosity, in our sample is 
star formation. 

The standard approach to estimating the dust atten- 
uation in star-forming galaxies is to use the flux ratios of 
Balmer lines. Under typical conditions (i.e. T c ~ 10 4 K), 
and in the absence of dust, the line ratio /h q //h^ = 2.86 
(Osterbrock 1989). Observed differences in this ratio, also 
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Figure 3. Composite rest-frame spectrum for the 57 24 n 250 fim 
sources with robust line detections in the region of Ho (top) and 
H/3 (bottom). The dashed line is the composite spectrum for only 
those 24 sources which are visible at the rest frame wavelength of 
H/3. The grey solid line represents the lcr variance in the compos- 
ite. Well-known spectral features are marked. Encouragingly, our 
composite spectrum recovers several other spectral lines, such as 
[Nil] 6584, the blended doublet of [SII] 6716 and [SII]6731, and 
weak signatures of [OIII] 5008 and H/3. 



known as the Balmer decrement, can be attributed to dif- 
ferential dust attenuation at the rest-frame wavelengths of 
the Balmer lines. 

Using our composite spectrum we estimate the aggre- 
gate value of the Balmer decrement for SPIRE sources at 
z Z 1. The H/3 line is visible in the FMOS wavelength 
coverage at 1.26 < z < 2.5, excluding sky lines. To ac- 
count for this we build a second composite spectrum, us- 
ing only those sources which have 'clean' (i.e. no overlap- 
ping sky lines) FMOS wavelength coverage at the wave- 
length of both H/3 and Ha; only 24 sources satisfy this cri- 
terion. The composite spectrum from these sources in the 
region of Ha and H/3 is shown in Fig. [3] The H/3 observ- 
able sources have a mean redshift of (z) — 1.36 and a mean 
IR luminosity of (Lib.) = 10 45 ' 5 ergss -1 (10 12 L Q ). We mea- 
sure line fluxes of (/Ha}=6.8±0.6 x 10 -16 ergscm _2 s _1 and 
(/h,s) = 1.3±0.4 x 10 -16 ergs cm~ 2 s~ 1 , after correcting both 
lines for stellar absorption (EWhq = 4.4 and EWh/3=2.8A; 
Moustakas & Kennicutt 2006) and applying an aperture 
correction of 2.8 (the mean value for these 24 sources; see 
i|3.3[) . This gives a Balmer decrement of (R) = /hq//h^ = 
5.2 ± 1.6, resulting in a dust attenuation of {E(B — V)) = 
log lo (i?/2.86)/0.4[fc(A Ha ) - kiXnfi)] = 0.51 ± 0.26, where 
fc(A Ha ) = 4.596 and fc(A Hj s) = 3.325 (Calzetti et al. 2000). 
This is equivalent to A v —2.1 mags, similar to that found 
for local IR-luminous galaxies (e.g. Hopkins et al. 2001, Wi- 
jesinghe et al. 2011). 



3.3 Relationship between Ha and IR star 
formation rate estimates 

In the absence of AGN activity or strongly non-solar metal- 
licity, differences between the Ha estimated SFR (SFRhq), 





1.5 


X 
CC 
LL 


1.0 


CO 

o 


0.5 


of 

LL 




GO 

o 


0.0 


a> 




o 


-0.5 




SFR„.SFR, R +SFR Hl 




1.0 



0.5 



> 
i 

m 

0.0 LU 



-0.5 



0.0 0.2 0.4 0.6 0.81.0 

E(B-V) BD 



0.0 0.2 0.4 0.6 0.81.0 

E(B-V) BD 



Figure 4. Left panel: Comparison of dust attenuation estimates 
from the Balmer decrement (E(B — V)bd) and from the ratio 
of SFRtot /SFRh q (E(B - V) m ) for SDSS galaxies. The dashed 
line represents E(B — V)ir = E(B — V)bd- The two estimates 
are excellent agreement, with a RMS difference of 0.17. Right 
panel: Comparison of dust attenuation estimates from the Balmer 
decrement (E(B — V)bd) ana from the ratio of SFRK09/SFRjj a , 
where SFRko9 is the combined Ha+IR SFR estimator presented 
in Kennicutt et al. 2009 (K09). Again, the dashed line represents 
E(B — V) m = E(B - V)bd It can be seen that the K09 SFR cali- 
bration gives much poorer agreement with the Balmer decrement 
estimated attenuations (RMS= 0.27). 



and the best estimate of the total SFR (SFR to t) can be 
attributed to the effect of dust attenuation. Thus the ra- 
tio SFRtot /SFRhq can be used as an estimator of the level 
of dust attenuation (e.g. Hopkins et al. 2001; Kewley et al. 
2002) . In order to calculate SFR to t and SFRh q the following 
steps were taken. 

IR luminosities are calculated by fitting template SEDs 
to the 24 /im and SPIRE data, then integrating the best 
fit template in the range 8-1000 fim. IR template SEDs are 
taken from Rieke et al. (2009). Both IR and Ha luminosi- 
ties are converted into SFR via the relations presented in 
Kennicutt (1998), assuming a Chabrier (2003) IMF, i.e; 



SFRir = 2.61 x 10 _44 Lir (ergs" 
and 

SFRh q = 4.61 x 10" 42 L Ha (ergs 



(3) 



(4) 



Finally, SFR to t is calculated by adding together the IR and 
Ha star formation rates, i.e. 



SFRt, 



SFRir + SFRhc 



While other combined Ha+IR estimators of SFR ex- 
ist (e.g. Kennicutt et al. 2009), these may not give good 
results for the class of IR-luminous galaxies (nor Herschel 
derived IR luminosities). We can confirm our assumption 
that SFRtot /SFRhq is a good tracer of dust attenuation by 
comparing to estimates of the attenuation from the Balmer 
decrement for our SDSS sample. All 156 of our SDSS com- 
parison sample have reliable (> 3er) H/3 line flux estimates. 
We produce estimates of the dust attenuation independant 
of the far-IR measurements using the Balmer decrement. 
Fig. [4] compares the dust attenuation estimated from the 
Balmer decrement {E(B — V)bd) to the infered value as- 
suming SFRtot /SFRhq is a good estimate of Ana- 

Fig. [5] shows the relationship between SFRhq and 
SFRtot for our sample of 57 sources detected at 24 /im, 
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Figure 5. Comparison of Ha to total (Ha + Far-IR) SFRs, 
assuming the conversions of Kennicutt (1998) with a Chabrier 
(2003) IMF. FMOS detected sources (red triangles) and the low- 
2 sample of 160 /im-selected SDSS galaxies (blue squares) are 
shown. The dashed line shows the best-fit log-linear relation to 
the data. 



250 ;Um and Ha. Fitting this observed correlation with a log- 
linear function results in the best fit; 



log 10 SFRhc = (0.82 ± 0.08) log 10 SFRir + 0.47 ± 0.11 (M 

Using the Calzetti et al. 2000 model for the variation 
of dust attenuation with wavelength we can use this result 
to convert the ratio SFRtot/SFrW to E(B - V). Taking 
our best fit correlation between these values we derive the 
following relationship between SFR and dust attenuation; 

E(B-V)m. = (0.135±0.06) log 10 SFR tot +0.35±0.08 (M yr" 

Fig. [6] presents dust uncorrected SFR to t/SFRH a vs. 
SFRtot- We also show in Fig. [S] our best fit relationship be- 
tween E(B — V) and SFR to t, and the empirical relationship 
for low-z IRA S galaxies as derived by Hopkins et al. (2001). 
These two relations show reasonable agreement, with the 
Hopkins et al. (2001) relation slightly below our simple log- 
linear fit. 

Alternative estimates of the mean (and standard devia- 
tion) of dust attenuation measured from the Balmer decre- 
ment of individual SDSS spectra, and the aggregate FMOS 
value from our composite spectrum f i|3.2[) are also shown 
in Fig. [5] Encouragingly the direct estimates of the mean 
dust attenuation are in good agreement with that inferred 
from SFRir/SFRhq, and both the dust attenuation-SFR 
relations. 

In Fig. [6] are also shown unlensed 850 /im selected 
sources (SMGs) at z Z 2, with Ha line measurements from 
Swinbank et al. (2004). Here we make use of the SFR(Ha) 
and L(FIR) quantities given in Table 2 of Swinbank et al. 
(2004), converting SFR(Hq) to the Chabrier (2003) IMF 
used here and calculating SFRir from L(FIR) using the 
equation given above. SFR(Ha) as quoted by Swinbank et 
al. (2004) includes corrections for slit-loss and so should 
be compatible with the values we derive from FMOS and 
SDSS data. Interestingly the z z 2 SMGs appear slightly 
above both the Hopkins et al. (2001) and our best-fit dust 
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Figure 6. Ratio of IR to dust uncorrected Ha-based SFR vs. IR- 
based star formation rate (SFR). The equivalent E(B — V) is also 
given, assuming, a Calzetti et al. (2000) attenuation curve. FMOS 
detected sources (red triangles) and the low-z sample of 160 (im- 
selected SDSS galaxies (blue squares) are shown. The dashed line 
shows the local relation between SFR and E(B — V) from Hop- 
kins et al. (2001). The E(B - V) determined from the Balmer 
decrement of our FMOS composite spectrum (Fig. [3}, and mean 
value for SDSS spectra are shown as cyan and orange circles, re- 
spectively. Also shown are the positions of z ~ 2 SMGs with Ha 
line measurements from Swinbank et al. (2004). Dotted lines show 
)the effect of the Ha detection limit (5cr) at the mean redshift for 
sources with that SFR to t for the SDSS (blue) and FMOS (red) 
samples. 



attenuation - SFR relation, suggesting that SMGs experi- 
ence enhanced attenuation. However it is worth noting that 
lwgsthe L(FIR) estimates for the SMGs come from pre-Herschel 
submm-radio estimates and hence may be overestimated 
(see Magnelli et al. 2012). Future studies with FMOS in 
this SFR range, as well as a re-assessment of the SMG pop- 
ulation with Herschel photometry will allow this trend to be 
confirmed. 

The modest Ha detection limits achievable with FMOS 
mean we will not recover IR-sources which have very large 
SFRtot /SFRh q - To quantify this we calculate the typical 
maximum observable limit of SFR to t/SFRH a as a function 
of SFRtot for both the SDSS and FMOS samples. For the 
SDSS sample we assume a detection limit of fn a = 5.8 x 
10 16 erg cm -2 s _1 , while for the FMOS sample we assume a 
detection limit of fHa = 1 X fO 16 erg cm -2 s _1 . In both cases 
these limits include a correction for the mean loss due to 
the limited aperture of the fibres (1.2 arcsec diameter for 
FMOS, 3 arcsec diameter for SDSS). While the maximum 
limit for the SDSS sample is significantly higher (~ 0.2 dex) 
than the observed values of SFRtot/SFRHa, the limits for 
the FMOS dataset appear quite close to the observed data 
points. 

Given ~ 65 per cent of our parent sample is undetected 
in Ha a potential explanation for this large incompleteness 
is a significant population of sources with SFR to t/SFRH Q 
above these selection limits. Hence the observed consistency 
with the Hopkins et al. (2001), and our Eqn. [S] may be a 
result of a bias towards low SFR t ot/SFRH Q . 

While we cannot rule out the existence of large 
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SFRtot/SFRHa sources (as we cannot detect them), we can 
estimate the observed completeness for our parent sample 
assuming our best-fit to the E(B — V)-SFRtot relation is a 
good description for the whole population. For each source 
in our parent sample of 168 24 PI 250 /im sources we first es- 
timate Lir, assuming the photo-z from Ilbert et al. (2009) 
and the SED fitting process described above. The Ha line 
flux is then predicted from Lir (assuming SFRir =SFR to t) 
using the best fit log-linear relation from Fig. [5] (with 0.5 
dex of intrinsic scatter) and the mean loss due to the fibre 
aperture (2.8). Applying the completeness curves from Ta- 
ble[5]we would expect a completeness of 35 per cent, in good 
agreement with the observed completeness of 34±4 per cent, 
and consistent with a similar assessment of the completeness 
presented in q3.ll 

If a large fraction of our parent IR-selected sample 
had dust attenuation levels (i.e. SFR to t/SFRHa) signifi- 
cantly above that predicted by Eqn. [5] we would expect 
much lower completeness in our FMOS observations than 
achieved. Given the good agreement between the observed 
and expected completeness we conclude that Eqn. [S] must 
hold for the bulk of IR- luminous sources at z ~ 1. 

3.4 The mass— metallicity relation for IR-galaxies 

A significant fraction (28/57) of our sample have robust 
(SNR> 3) measurements of the [Nil] 6584 line. This allows 
us to investigate the gas phase metallicity (12+log 10 O/H) 
of our IR-sources via the ratio of the [Nil] 6584 to the Ha 
line (N2 method; Kewley & Dopita 2002; Pettini & Pagel 
2004). For each source with a robust [Nil] measurement we 
estimate the metallicity via Eqn. 1 of Pettini & Pagel (2004; 
PP04); 

12 + log 10 O/H = 8.9 + 0.57 x log 10 / [N n]// H a 

The existence of a relationship between stellar mass and 
metallicity has been confirmed across a wide range in red- 
shift (0 < z < 3; Lequeux et al. 1979; Tremonti et al. 2004; 
Erb et al. 2006; Maiolino et al. 2008; Zahid et al. 2011), with 
a steady trend towards lower metallicity for a given stel- 
lar mass with increasing redshift. To investigate the mass- 
metallicity relation for our z ~ 1 IR-sources we combine 
our metallicity estimates with stellar masses as derived by 
Wang et al. (2012). Stellar mass estimates from Wang et 
al. (2012) are calculated by finding the best-fit stellar popu- 
lation model to the observed multi-wavelength photometry 
using the Le Phare software (Arnouts et al. 2002; Ilbert et 
al. 2006), combined with stellar population synthesis mod- 
els from Bruzual & Chariot (2003) and assuming a Chabrier 
(2003) IMF. 

The left panel of Fig. compares the stellar mass to the 
metallicity for both our sample of z ~ 1 Herschel-FMOS and 
comparison sample of z ~ 0.1 SDSS sources. A tentative 
trend of metallicity with stellar mass can be seen in both 
samples, although with significant scatter. No discernable 
evolution is seen in the metallicity between the z ~ 1 and 
z ~ 0.1 IR-selected samples. 

Also shown in the left panel of Fig. [7] are studies using 
stacking of near-IR spectra by Erb et al. (2006) with Keck 
NIRSPEC atz~ 2.2 and Yabe et al. (2012) with FMOS 
at z ~ 1.4 (We convert the Yabe et al. 2012 stellar masses 
to a Chabrier 2003 IMF for consistency). Given the similar 



redshift range between our Herschel-FMOS sample and the 
Yabe et al. (2012) study the difference in the metallcities is 
surprising, although likely due to a combination of the dif- 
ferent mass ranges probed, the stacking nature of the Yabe 
et al. (2012) result and our bias towards dusty galaxies. 

In the middle panel of Fig. [7] we show the relation- 
ship between specific SFR (SSFR;SFR/M*) and metallic- 
ity. Over the range in SSFR well-sampled by our Herschel- 
FMOS and SDSS samples we see no discernable trends, al- 
though there are hints of decreasing metallicity with increas- 
ing SSFR above SSFR= 0.5 Gyr~ x . It has been proposed 
that the metallicity is a function of both SFR and stellar 
mass, with a "fundamental mass relation" (FMR) linking 
the three parameters (Mannucci et al. 2010). In Fig. [7] we 
show the fit to the FMR from Mannucci et al. (2010) for 
values of SFR representative of our study. The metallcities 
of both our Herschel-FMOS and SDSS are consistent with 
the prediction from the FMR. 

Finally, we consider the relationship between dust at- 
tenuation and metallicity. In the right panel of Fig. [7] is 
shown the metallicity as a function of dust attenuation 
(as traced by the ratio SFRtot/SFRHa)- A clear trend be- 
tween dust attenuation and metallicity is seen in both the 
//er,sc/ie/-FMOS and SDSS samples, with the most metal- 
rich galaxies experiencing the largest levels of attenuation. 
The best fit log-linear relation between attenuation and 
metallicity is found to be: 

12 + log 10 O/H = 0.19 log 10 SFRtot /SFRh q + 8.57. 

Converting to E(B - V) via the Calzetti et al. (2000) 
attenuation curve we find Eqn. [6] 

12 + log 10 O/H = 0.24 E(B-V) + 8.57 (6) 

This correlation is somewhat expected; the dust grains 
responsible for attenuating starlight (and nebular line emis- 
sion) are synthesised from metals in the ISM. Similar stud- 
ies, using the IR to UV luminosity ratios as a proxy for dust 
attentuation, have also found a correlation between metal- 
licity and attenuation for nearby galaxies (Heckman et al. 
1998) and z ~ 2 (Reddy et al. 2010). Given the correlation 
between metallicity with attenuation the high metallicity of 
our Herschel-FMOS sample (when compared to Yabe et al. 
2012 at z ~ 1.4) is to be expected as our IR-selected sample 
must be biased towards the most obscured sources. 



4 DISCUSSION AND CONCLUSIONS 

We have investigated the properties of z ~ 1 IR luminous 
galaxies by performing near-IR spectroscopy of a sample of 
Spitzer and Herschel selected sources with FMOS. Candi- 
date emission lines were identified in the 2D reduced FMOS 
frames via a semi-automated procedure. Via comparison 
with known spectroscopic redshifts, and direct testing of the 
line detection algorithm, we estimate that our Ha line sam- 
ple is z 90 per cent reliable. Our scientific conclusions are; 

• Robust detections of Ha were found for 57 of 168 
24P1250 /im-sources, resulting in a detection rate of 34±4 per 
cent. This detection rate is consistent with the expected in- 
completeness measured by simulating the line detection pro- 
cess. For sparse targets such are i/ersc/ieZ-selected sources, 
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shown for those sources without robust [Nil] line measures. Sources from our Herschel-FWlOS sample at z ~ 1.2 are shown as red dots, 
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& Ellison (2008); KE08) and z ~ 2.2 (Erb et al. 2006; E06) are shown as dotted and dot-dashed lines in the left panel, respectively. 
Results from a previous FMOS study for star-forming galaxies at z ~ 1.4 (Yabe et al. 2012) are shown as orange stars. Where necessary, 
comparison samples have been converted to use the Pettini & Pagel (2004) calibration of the N2 metallicity tracer (KE08) and a Chabrier 
(2003) IMF. In the middle panel the position of the "fundamental metallicity relation" of Mannucci et al. (2010) is shown for a range 
of SFRs representative of our sample as the dotted lines. In the right panel the best-fit log-linear relation between metallicity and dust 
attenuation (as traced by the ratio SFRtot/SFRna) is shown as a dashed line. 



FMOS is a competitive redshift recovery instrument with 
equivalent optical MOS instruments on other 8m class tele- 
scope (e.g. VIMOS on VLT). 

• The mean dust attenuation, estimated via the Balmer 
decrement for the Ha and H/3 emission lines for a composite 
spectrum of Ha detected sources, is E(B — V) = 0.51 ± 0.27 
for Lir = 10 12 Lq sources at (z) = 1.36. 

• Good agreement was found between the dust attenu- 
ation estimated from the Balmer decrement and that in- 
ferred from the ratio of Ha-estimated to best estimate of 
the total star formation rate. Using SFRtot/SFRHa as an 
indicator of attenuation in the Ha line we derive a relation- 
ship between dust attenuation and star formation rate of 
E(B -V) = (0.135 ± 0.06) log 10 SFR tot + 0.35 ± 0.08. These 
results are broadly consistent with the relationship between 
star formation rate and dust attenuation seen both in \aw-z 
star forming galaxies (e.g. Hopkins et al. 2001). 

• The gas phase metallicity relation was investigated for 
the subset of Herschel-FMOS, and SDSS, sources with ro- 
bust measurements of the [Nil] 6584 line. For IR-selected 
sources with M* ~ 10 10,5 Mq the typical metallicity is found 
not to evolve between z ~ 0.1 to z ~ 1.2. No discern- 
able trend with specfic SFR is seen, but a strong correla- 
tion between metallicity and dust attentuation is seen, de- 
scribed by a best-fit log-linear relation; 12 + log 10 O/H = 
0.24 E(B - V) + 8.57. 
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Tabic Al: Details of Herschel 250 /im and Spitzer 24 /^m-selected sources with reliable Hct 
line detections. 
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Tabic Al: (continued) 
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15.70 ±2.74 


< 12.59 


< 17.81 


4.58±0.46 


1.10±0.31 


150.31065 


2.31922 


1.43 


25.53±0.13 


0.17±0.01 


< 8.21 


< 18.13 


< 19.54 


3.28±0.26 


< 2.85 


150.37685 


2.45953 


1.30 


25.53±0.14 


0.36±0.01 


19.29 ±2.74 


22.76 ±4.10 


< 19.70 


4.51±0.38 


< 2.60 


150.46438 


2.63802 


1.46 


24.59±0.10 


0.09±0.01 


17.27 ±2.74 


25.30 ±3.15 


26.53 ±4.91 


1.80±0.32 


< 1.50 


150.68825 


2.45364 


1.38 


25.95±0.15 


0.18±0.01 


19.43 ±2.73 


16.53 ±3.66 


< 17.36 


6.78±0.45 


2.20±0.27 



a Position from -Zf814W catalogue. 

b Redshift from detected Ha: line. 

c Limits are 3<r, including confusion noise. 



d All values corrected for the limited aperture (1.2 arcsec diameter) of the FMOS fibre. 
e All values corrected for the limited aperture (1.2 arcsec diameter) of the FMOS fibre. Limits arc 3ct, no value given in cases where the [Nil] 6584 line falls on a OH sky line. 



